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Abstract. 

" Q I Recently hybridized monolayers consisting of hexagonal boron nitride (ft,-BN) 

CZJ ■ phases inside graphene layer have been synthesized and shown to be an effective way 

of opening band gap in graphene monolayers [1]. In this paper, we report a first- 
principles density functional theory study of the /i-BN domain size effect on the elastic 
properties of graphene/boron nitride hybrid monolayers (/i-BNC). We found that both 
in-plane stiffness and longitudinal sound velocity of /i-BNC linearly decrease with h- 
BN concentration. Our results could be used for the design of future graphene-based 
nanodevices of the surface acoustic wave sensors and waveguides. 
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1. Introduction 

Substituting C with B and N atoms in graphene has been shown to be a promising way 
to improve semiconducting properties of the graphene [2l |3l H] . Hexagonal boron nitride 
(^-BN) monolayer [5] and graphene have similar 2D lattice structures but with very 
different physical properties. Interesting nanostructures can be made by mixing these 
two structures[6l [3 El E]. 

Recently, a promising method has been reported in which atomistic monolayers 
have been generated consisting of /i-BN phases in graphene (/i-BNC) using a thermal 
catalytic chemical vapor deposition method [Tj. These hybrid monolayers have been 
shown to have isotropic physical properties which can be tailored by controlling the 
kinetic factors affecting the /i-BN domain size within graphene layer. This is different 
from B-doped or N-doped graphene, where the integrity of the /i-BN structure is missing. 

The electronic band structures of the /i-BNC heterostructures were studied in 
our previous work pj]J. However, the mechanical properties of these heterogeneous 
nanosheets are still unknown. In this paper we investigate the elastic properties of the 
/t-BNC monolayers as a function of /i-BN concentration using ab initio density functional 
theory. 

2. Modeling and Computational Details 

The effect of /i-BN domain on the properties of the /i-BNC hybrid structures is modeled 
by only considering the eect of its size while maintaining its hexagonal structure within 
the system. The h-BN domain size effect then can be represented by the h-BN 
concentration x in the model as (B3N3)j.(C6)i-a; where (B3N3) and (Cg) denote the 
nanodomain structure of h-BN monolayer and graphene, respectively. The proposed 
domain size effect model (B3N3)a,(C6)i-x is based on the results of previous studies 
of Bj;CyNz layered structures [H], layers and nanotubes[T2l [13], quantum dots and 
nanorods[14j. and monolayer nanohybrids[15j . It is a general belief that the /i-BN 
segregates in the /i-BNC, and the system gains lower energy, larger band gap, and 
better thermal stability after phase segratation [H [121 [IS]- The six-atom hexagonal 
structures (both B3N3 and Cg) are the basic blocks in these heterostructures. This 
model of (B3N3)2,(C6)i-x captured the main feature of these heterostructures. As shown 
by previous works pT| [12] . both stoichiometry and geometry change the band gap of 
B-C-N nanotube. The maximum band gap is achieved at B/N ratio of 1. In the B-C- 
N monolayer, however, the domain size is a dominant factor comparing with domain 
shape in changing the band gap [15]. By varying h-BN concentration x, the domain size 
effect on mechanical properties of mono-layer hexagonal BNC heterostructures can be 
appropriately presented in our (B3N3)^.(C6)i_a. model. 

We need to emphasis that although we only investigated the effect of stoichiometry 
of /i-BN here, this domain model is different from point model of Ba-C^N^ model at 
X = z where hexagonal (B3N3) structure is not considered [H]. In other words, our 
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Figure 1. Atomic structures of five configurations in order of /i-BN concentration: 

(a) 0%; (b) 25%; (c) 50%; (d) 75%; (e) 100%. 



model specified not only the stoichionietry of B/N ratio of 1, but also the hexagonal 
(B3N3) and Cq structures, to represent the two separated phases in heterogeneous h- 
BNC structures. 

We examined the change in elastic properties of the /i-BNC monolayer as a function 
of /i-BN concentration. Five /i-BNC configurations, in order of /i-BN concentration, 
0%, 25%, 50%, 75% and 100% have been studied, where 0% and 100% corresponding 
to pure graphene and /i-BN, respectively. The three other concentrations were selected 
based on their simplicity and representativeness. The atomic structures of these five 
configurations (Fig. [1]) were determined by the ab initio density functional theory 
through geometry optimization. 

DFT calculations were carried out with the Vienna Ab-initio Simulation Package 
(VASP) [ini [171 [13 [19] which is based on the Kohn-Sham Density Functional Theory 
(KS-DFT) [201 [21] with the generalized gradient approximations as parameterized by 
Perdew, Burke and Ernzerhof (PBE) for exchange-correlation functions [22]. The 
electrons that explicitly included in the calculations, are the (2s^2p^) electrons of carbon, 
the (2s^2p^) electrons of boron and the {2s'^2p^) electrons of nitrogen. The core electrons 
(Is^) of carbon, boron and nitrogen were replaced by the projector augmented wave 
(PAW) and pseudopotential approach|23l |2l]. A plane-wave cutoff of 520 eV was used 
in all the calculations, such value was chosen to eliminate the pulay stress during the 
geometry optimization [25]. 

The criterion for stopping the relaxation of the electronic degrees of freedom was 
set by total energy change to be smaller than 1.0—6 eV. The optimized atomic geometry 
was achieved through minimizing Hellmann-Feynman forces acting on each atom until 
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the maximum forces on the ions were smaller than 0.001 eV/A. 

The atomic structures of the five configurations were obtained by fully relaxing a 
24-atom-unit cell where all atoms were placed in one plane. The irreducible Brillouin 
Zone was sampled with a Gamma-centered 19 x 19 x 1 /c-mesh and initial charge densities 
were taken as a superposition of atomic charge densities. There was a 14 A thick vacuum 
region to reduce the inter-layer interaction to model the single layer system. 

The elastic tensors were determined by performing six finite distortions of the lattice 
and deriving the elastic constants from the stress-strain relationships [26] . In this study, 
we were interested only the in-plane elasticity. The strains of ±0.015 were applied along 
two directions, x and xy. Two stress-strain curves, an vs en and ai2 vs £12, were 
computed, where an and Sn are stress and strain in xx direction, respectively, and ai2 
and £12 are the stress and strain in the xy direction. The elastic constants then were 
obtained through the least-squares extraction of coefficients from calculated stress-strain 
data. The final results of the values of these elastic constants were printed directly in 
VASP [26]. 



3. Results and Analysis 

With the full geometry optimization, including allowing the change of shape, the final 
atomic structures of the five configurations were determined using DFT calculations, as 
shown in panels (a)-(e) in Fig. [H The hexagonal structures were still kept, even for the 
configuration (c), which has 50% /i-BN. In addition, we tested another atomic structure 
of configuration (c), where the two gaphene rings are in diagonal, in stead of parallel 
along the side. The lattice constants of the two atomic structures are the same after 
geometry optimization, since the two atomic structures are symmetrical, and both have 
one layer of graphene ribbon in addition to one layer of /i-BN ribbon, due to the periodic 
boundary conditions. 

The lattice constants of the /i-BNC mixtures were calculated as the half of the lattice 
vectors of the super-cells, for the convenience of comparison to those of the pure graphene 
and /i-BN. We found that the lattice constant increases with /i-BN concentration, which 
is denoted as Cbn in this paper. Our results are in a good agreement with previous 
experimental results reported for /i-BN (2.51 A) [27] and graphene (2.46 A) [28] . 

The elastic constants were obtained from DFT calculations. Due to the symmetry, 
only Cii and C12 are independent. Cn decreases linearly with respect to Cbn^ as 
plotted in top-left panel of Fig. O The in-plane stiffness Ys can be obtained from the 
elastic moduli Cn and C12 as Yg = {Cfi — Cfg)/^^!!- The Poisson's ratio z/ which is the 
ratio of the transverse strain to the axial strain can be obtained from elastic moduli as 
u = C12/C11. Our result of u and C for the five configurations are shown in bottom 
panels of Fig. |2l Our calculated value for in-plane stiffness of graphene (347.2 N/m) is 
in a good agreement with the experimental value (340 ± 50 N/m) [29], and theoretical 
predictions (348 N/m in ref.[30] and 335 N/m in ref. [31]). Our calculated value of Yg 
for /i-BN (291.3 N/m) agrees with ah initio (GGA-PW91) prediction (267 N/m in ref. 
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Figure 2. Elastic constants Cn and C12, Poisson's ratio v and in-plane stiffness Ys 
as a function of /i-BN concentrations C'bn- 



|31]). The calculated Poisson's ratio u is 0.16 for graphene and 0.18 for /i-BN, which 
agree with reported values of 0.16 and 0.21, respectively [3T] . 

As shown in Fig. |21 both the Cu and the in-plane stiffness decrease nearly linearly 
as /i-BN concentration increases. However, the C12 and Poisson's ratio, show rather 
complicated behaviors. The similar trend of the Cu and the in-plane stiffness is due to 
the fact that Cu is the dominant factor in computing the in-plane stiffness, about six 
times bigger than C12. For the same reason, Poisson's ratio u has the similar trend as 
C12. Therefore, our analysises only focus on the dominant factor Cu. 

The mechanical behaviors are ultimately determined by the electronic structures. 
For a better understanding of the mechanism of linear elastic properties with respect to 
the Cbn, we studied the total electronic charge density pe = q/S, where q is the total 
electronic charge of the system, and S the area (cross section) of the system in the basal 
plane. We found that pe monotonically decreases with Cbn, as plotted in the top-right 
panel in Fig. [31 This relationship is consistent with the q-Cbni since q is the same in 
all five configurations. The area mass density pm = m/S (m is the atomic mass) are 
also studied, as plotted in the bottom-left panel in Fig. |3l pm reaches the minimum at 
Cbn = 0.5. 
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Figure 3. Lattice constants a, electronic cliarge density p^ and mass density pm as a 
function of /i-BN concentrations Cbn- The in-plane stiffness Ys varied with electronic 
charge density is plotted. 



The relationship between Yg and pe was examined and plotted in bottom-right panel 
of Fig. 131 It is interesting to note that the in-plane stiffness monotonically increase with 
respect to the area charge density pe- It gives a hint that the elastic properties can be 
engineered by doping or introducing defects which changes the charge densities. 

In these /i-BNC structures, there is a non-zero stiffness both for volumetric and 
shear deformations. Hence, it is possible to generate sound waves with different velocities 
dependent on the deformation mode. Sound waves generating volumetric deformations 
(compressions) and shear deformations are called longitudinal waves (p— wave) and 
shear waves (s— wave), respectively. The sound velocities of these two type waves are 
respectively given by: 



(1) 



YJl 



p„(l + z/)(l-2z/) 



'G 



12 
Pm 



(2) 



The dependence of Vp and Vg on Cbn are plotted in Fig. HI Vp monotonically reduced 
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Figure 4. p— wave and s—wave velocity as a function of /i-BN concentrations Cbn- 



with Cbn- However, Vg has a more complex behavior. The Vg of /i-BNC decreases with 
Cbn^ but is lower than that of the pure phases, both graphene and /i-BN. The Vp/vg is 
in the range between 2.42 and 2.64. 

Since the speed of these sound waves can be measured experimentally when these 
/t-BNC heterostructures are synthesized, the predicted p— wave and s— wave velocities 
could be served as the quantities for the validation of the elastic properties. 

Furthermore, we predicted a linear relationship between the sound velocity and 
the concentration of /i-BN. As a result, a sound velocity gradient can be achieved by 
introducing /i-BN domains into graphene. The sound velocity gradient could be used 
to form a sound frequency and ranging channel, which is the functional mechanisms of 
surface acoustic wave (SAW) sensors and waveguides. Thus, graphene-based nano- 
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devices of SAW sensors and waveguides can by synthesised with /i-BNC for next 
generation electronics. 

4. Conclusion 

In summary, we used ab initio density functional theory to investigate the effect of 
the /i-BN domain size on the elastic properties of graphene//i-BN hybrid monolayers. 
The elastic constants of five configurations were explicitly examined. We found that 
the in-plane stiffness increases linearly with the /i-BN concentrations. We predicted 
a linear relationship between the sound velocity and the concentration of /i-BN. This 
knowledge could be used for the design of future graphene-based nanodevices of the 
surface acoustic wave sensors and waveguides. Thus our results may provide guidance 
in practical engineering applications of these nano-heterostructures. 
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